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Abstract 

An effective action for Bose-Hubbard model with two- and three-body 
on-site interaction in a square optical lattice is derived in the frame of a 
strong-coupling approach developed by Sengupta and Dupuis. Prom this 
effective action, superfluid-Mott insulator (MI) phase transition, excita- 
tion spectrum and momentum distribution for two phases are calculated 
by taking into account Gaussian fluctuation about the saddle-point ap- 
proximation. In addition the effects of three-body interaction are also 
discussed. 



PACS number(s): 03.75.Hh, 73.43.Nq, 05.30.Jp 



1 Introduction 

The quantum phase transition from superfluid to Mott insulator (MI) has been realized for 
cold bosonic atoms in a cubic optical lattice [Ij and well explained in the frame of Bose- 
Hubbard model (BHM) only with on-site two-body interaction [2]: when the lattice potential 
is very weak, the kinetic energy dominates and the system is a delocalized superfluid. By 
contrast if the lattice potential is very strong, the interaction will forbid bosons to hopping 
from one site to another site and the system is in the MI phase. Hence for intermediate 
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lattice potential we naturally expect the phase transition from superfluid to MI phase. In 
fact many other phenomena of bosons in optical lattices can also be well described by BHM 

However, the work of Kohler [5] and Johnson et al. [6j pointed out that an effective 
three-body interaction should be included to explain the experimental criterion for collapse 
and revival of coherent matter waves |8l [9]. In terms of a pure bosonic system, the 
effective three-body interaction is related to the two-body scattering length and the 
dilute gas parameter r/ = Vpaf by gs oc In [Cr;^] [SJ [lOl [TTl |T2], where p is the boson 
density and the constant C has been determined by applying a microscopic description [5]. 
Generally in the current experiments the strength of three-body interaction is much weaker 
than its two-body counterpart. In order to enhance the ratio of three-body to two-body 
interactions, a feasible route is by the induction of fermions. Considering the mixture of a 
BEC and a single component Fermi gas, the Fermi degrees of freedom can be completely 
integrated out and the effective three-body interaction can be produced [HI Ej. In this 
method the interspecies interaction gsF not only induces three-body interaction g^ oc g%p, 
but also weakens two-body interaction to g2 = gs — PQbf^ where gs and g2 are bosonic 
two-body interaction before and after the integration of Fermi degrees of freedom and P is 
a constant [13]. Very explicitly by Feshbach resonance to tune interspecies interaction gsF 
[13 dS], one should have enough space to adjust the 53/52 • From this viewpoint it is very 
significant to include three-body interaction into BHM. 

As a matter of fact, some references [T71 [181 UHl HO] have studied superfluid- MI phase 
transition in the presence of two- and three-body on-site interactions at the mean field level 
and found some interesting phenomena such as the rotation of phase boundary around a 
fixed point and extension of MI. Especially the reference [19] also studied the excitation 
spectrum in MI by functional-integral method. Different from these works, in this paper, we 
discuss the effects of three-body interaction on excitation spectrum and momentum distri- 
bution both in the superfluid and MI phases on the same footing, following a strong-coupling 
approach [211 122] . From the experimental viewpoint, the momentum distribution can be 
measured by imaging atom gas after turning off the trap potential, while the excitation 
spectrum resulting from a particle-hole excitation can be attained by applying a potential 
gradient to the system in the MI phase. It is by probing the momentum distribution and 
excitation spectrum that superfluid-MI phase transition is observed [Ij. In section 2 we 
derive an effective action by performing two successive Hubbard-Stratonovich transforma- 
tions of the intersite hopping term. In section 3 starting from this effective action, the 
phase diagram, excitation spectrum and momentum distribution are calculated by taking 
into account Gaussian fluctuations about the mean-field approximation as in the Bogoli- 
ubov theories of the weakly interacting Bose gas. A brief conclusion is given in section 4. 
Throughout this paper, we set h = ks = 1, and take the lattice constant unity. 
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2 Hamiltonian and the Effective Action 



The system of bosons in an optical lattice with on-site two- and three-body interaction can 
be described by the modified Bose-Hubbard model 

H = - Y^tijalaj + ^^ni(ni - 1) + -yX! ~ ^)("-« ~ 2) ~ ^X!"!"*' 

i,j i i i 

where Oj is the bosonic field operator, rii = a\ai is the particle number operator for bosons 
at the lattice site i. The intersite hopping matrix tij is nonzero (labeled t) only when lattice 
sites i and j are the nearest neighbors, U , W are two- and three-body repulsive interaction 
strength among bosons. In the light of the statement in the introduction about the strengths 
of two- and three-body interactions, one consider that the ratio W/U is adjustable in a large 
parameter space. The last term involving the chemical potential is added because it is 
very convenient in the grand canonical ensemble. Without loss of generality a square optical 
lattice is assumed. 

In the imaginary-time functional integral formalism, the partition function of the system 
is written as 1 23 1 

Sq = (ir^a*(c?r - /u)aj + ■^^ni(ni - 1) + ^ ^ni(ni - l)(ni - 2), 



S' = dT 







-t a|aj 



<i,j> 



(2) 



with /3 = 1/T and r imaginary time. &[■■■] represents the functional integral for field 
operators. 

Carrying out Hubbard-Stratonovich transformation for S' and integrating out a, the 
partition function 



Zo [ ^[6*,6]e-^'**«^'''^+^[^*'''l (3) 



where denotes the inverse matrix of the intersite hopping matrix tij, Zq = J ^[a*,a]e~^° 
is the partition function in the local limit (t = 0), < • • • >o means that the average is taken 
with Sq. According to the linked cluster expansion [23], the functional V[b*, b] only includes 
contributions from linked diagrams. Until the fourth order of b and after a careful calculation 

V[b*,b] = - dndT2j2b*iri)Glc{Tl,T2)biiT2) + 

Jo 
1 fi^ 

- dTidT2dT3dT4 ^ G2c (n , T2 , T3 , r4)6* (ti)6* (T2)6i (r4)6i (ts) , (4) 
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where Gic{ti,T2) (G2c(ti, r2, ts, r4)) denote the connected local single (two) particle(s) 
Green function. In fact, is the generating functional of connected Green functions 

in the local limit if b field is regarded as external source. 

It is easily found that the relation between the Green function of a field Ga and Gb of 
b field is Ga = — t^^G^t^^. Hence once Gi, is known, we can get Ga- Unfortunately 
Green function Ga obtained in this way is not physical since it leads in the superfluid phase 
to a spectral function which is not normalized to unity [21]. To overcome this difficulty, a 
second Hubbard-Stratonovich transformation was performed to decouple intersite hopping 
term and the partition function becomes 

?[c*, c, 6*, 6]e<*»^-^^-<^'-^''=»+^[''*'^l (5) 

Reference [21] has shown that the auxiliary field c has the same correlation functions as the 
original boson field a, hence below we will substitute a for c. In ([5]) considering the second 
order term in y[&*,6] as weight and integrating out b field, till the fourth order of a, the 
effective action 5 [a*, a] is 



rl3 

S[a*,a] = - dTidT2^a*{Ti) G{} {ti,T2) + tij6{Ti - T2) aj{T2) 

+- / c?riciT2(iT3dr4 Vr(ri,T2,T3,r4)a*(Ti)a*(T2)aj(T4)ai(T3), (6) 
4 Jo , 



where we have discarded the contributions from all anomalous terms in the same spirit of 
Ref . [22] . r(Ti, T2, r3, r4) is the exact two-particle vertex in the local limit 

r(ri,T2,T3,r4) = - dT[dT2dT3dT!^G2ciT[,T2,T3,T4) X 

Jo 

Gri(Ti,rOGri(r2,r^)Gr,^(T4,Tl)Gri(T3,r^), (7) 

and has a complicated time-dependence. Taking the static limit of T into consideration and 
introducing a parameter 

— Ir — 

9 — i^i' static — „ — 4 5 {") 

where Gic, G2C denote the zero-frequency component of Fourier transformation of Gic, G2C 
respectively, we finally obtain 

f/3 



5[a*,o] = - dTidT2^a*{Ti) G^} {ti,T2) + tijd^n - T2) aj(r2) 
° id 

■ / dr^^ a*a'^aiai. (9) 
Jo . 



2 Jo 



In contrast to the original action on the one hand the intersite hopping has been 
renormalized by the exact local single particle Green function Gic, on the other hand the 
interaction U has been substituted by the exact local two-particle vertex. This action ([9]) 
is the starting point of our analysis. 
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3 Excitation Spectrum and Momentum Distribution 



In this section we decide the boundary of superfluid-MI phase transition, excitation spec- 
trum and momentum distribution for both superfluid phase and MI at zero temperature. 
Before doing these, we first calculate local Green functions Gic and G2c at zero temperature 
to determine the parameter g. 

In the local limit, the Hamiltonian ([1]) is diagonal about lattice sites and its eigen- 
states are Fock states \n >= (n!)~-'^/^(a'^)"|0 > with eigenvalue e„ = (C//2)n(n — 1) + 
(W/6)n{n — l){n — 2) — fin. In order to attain g we firstly calculate one-particle Green 
function G'i(r) = — < 7^a(r)a"'"(0) > and two-particle Green function G2(ti, T2; rs, r4 = 
0) =< 7^a(ri)a(r2)a'f (0)a"^(r3) >, which are easily calculated using the closure relation 
Y.^=o\n >< n\ = 1. Then Gic(r) = Gi(r) and G2c(ri, rs; rg, = G2(ri, rs; rg, - 
Gi{ti, T3)Gi{t2, T4) — Gi{ti, T4^)Gi{t2, T3) . Making Fourier transformation, we obtain 

Gic{tOJn) = -. h -. , (10) 



G 



2c 



4(no + l)(no + 2) 4no(no - 1) 



4no(no + I) _^ 4no(no + 1) 



(^710 ^no+l){^no £no-l)^ (^no £710+1 )'^ (^710 £710-1) 

^ H I 4(no + l)^ ^^^^ 

(£710 ~ £710-1)"* (£no ~ £710+1)"* 

where no is the occupation number of ground state in the local limit which minimizes the 
eigenvalue e^- 

3.1 Superfluid-MI Phase Transition 

To proceed further, we consider a quadratic expansion of the action ([9|) in terms of fluctu- 
ation near the saddle-point value of a field. Choosing saddle-point value Vo for that is 
o-iiT) — ^ V'o + c^iiT), dH becomes S = Sq + Si + S2 with 

50 = -(GIc+^t)lPl + ^^Pi 

51 = -MGic + 4t - gi^l) f dr ^[ai(r) + a*(r)], 

Jo i 

52 = - I (iridr2 ^a*(ri) Gi}{ti,T2) + tij5{Ti - T2) aj{T2) 

g rl3 

/ dry^ a*a* + aiUi + 4a* a^. (12) 
2 Jo Y 

The constant ipQ is determined by requiring the coefficients of the linear terms in Sj and a* 
to vanish, leading to the equation 

Gi,^ +At-g^Pl = 0. (13) 
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Figure 1: The phase diagram of superfluid-MI phase transition with different three-body 
interaction W. 



It is weh-known that a nonzero ipQ signifies superfluid of system, therefore the boundary of 
superfluid-MI phase transition is 



+ 4t = 0, 



(14) 



which is in agreement with the results in \n\ [T8| \T9\ I20j . Fig. 1, the phase diagram in the 
t/U—(j,/U plane for different W/U from (jl4p . shows that the effects of three-body interaction 
on the phase boundary are very dramatic. On the one hand Mott lobes with no = 1 remains 
unaltered in the presence of W . On the other hand for Mott lobes with higher densities, 
with the increase of W , their widths and heights gradually increase. In addition for small 
W phase boundary has the same characteristics as the conventional superfluid- MI phase 
transition that Mott lobes with higher occupation have smaller area. While W gradually 
increases and dominates two-body interaction [/, Mott lobes with higher occupation have 
bigger area. These phenomena can be illustrated from the pure two-body and three-body 
interaction |17) . 

The excitation spectrum and momentum distribution can be deduced from the second 
order action 5*2. Introducing the Fourier transformation 



a{k,iujn) 




(ir^ai(r)e-*^-^»+^'^"^, 



(15) 



where is the total number of lattice sites, Ri is a lattice vector for site i, Un = 2mT/f3 is 
bosonic Matsubara frequency, S2 becomes 



So 



a*{k,iujn),ai-k,-iujn) 



with 



M- 



2x2 



9^0 



M2x2 a{k,iujn),a*{-k,-iu}n) 



(16) 



(17) 
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(b) k 



Figure 2: The excitation spectrum for quasiholes (a) and quasiparticles (b). We choose 
t/U = 0.01, /i/C/ = 1.5 and W/U = 0,2.0,4.0 to make the system into MI with the filhng 
factor TiQ = 2. In (a) the effects of different W is very small. In (b) from bottom to top W 
is increasing. 

and Ek = —2t{coskx + cos/cy) denoting the boson dispersion in the absence of on-site 
interaction. 



3.2 MI Phase 

In the MI phase V'o = 0) tl^6 matrix M2x2 becomes diagonal and 5*2 is reduced into 
S2 = ^a*ik,iujn) -G^j{iujn) + £k a{k,iu}n)- 

k,UJn 



(18) 



The Green function gMiik,iu}n) = — < a{k,iujn)a*ik,iujn) > can be directly obtained 
gllj{k,iu!n) = G^^{iujn) - £k- Using ([10]), we obtain 



9Mi{k,iuj„ 



Zk 



+ 



I- Zk 



iun - El iujn - 



with 



Et 
Ak 



Zk 



+ ek±VAk 



(eno+l — Eno-l + £fc) — 4(en,)_i — eno)(eno — £no+l) 

no ^no+l ) + (no + l)(e„o - eno-i)], 
E^ - no{eno - enp+i) - {np + l){eno ~ ^no-i) 

K - El 



(19) 

(20) 

(21) 
(22) 



In (j20p . El stand for excitations of quasiholes (removing particles from sites) or quasi- 
particles (adding particles to sites), which is consistent with the results in [T7] and is shown 
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3-3 ' S 

Figure 3: The excitation spectrum for density fluctuation H^. We choose t/U = 0.06, 
W/U = 2.0 and /x/C/ = 1.5,5.5, 10.5 to make the system into MI with filling factor no = 
2, 3, 4 respectively. From bottom to top ^ is increasing. 



in Fig. 2(a) (Fig. 2(b)). Very explicitly, three-body interaction W has little impacts on quasi- 
hole excitations, but it has the large effects on quasiparticle excitation: with the increase 
of W , quasiparticle energy becomes higher, hence its excitation also becomes difficult. 

In Greiner's experiment [T], Superfluid-MI phase transition is also detected by applying 
a potential gradient to the system in the MI phase and probing the excitation spectrum 
resulting from a particle-hole excitation. According to the above results, we can find a 
first approximation for the dispersion of particle-hole excitation (density fluctuations) 
by subtracting E'^ from E'^ , which yields 

Hfc = E+ - E^ = VA;. (23) 

Due to little impacts of three-body interaction W on £'^, the dependence of on W is the 
same as E^ . In Fig. 3 we show for a fixed W , but different filling factor tiq, i.e. different 
chemical potential ^. We can see that in MI phases there is always a band gap increasing 
with the filling factor, which is different from the situation without three-body interaction 
where the gap is approximately independent of the filling factor. It is from this perspective 
that reference |20] proposed to observe weak three-body interaction for high filling MI. 

From the theorem of fluctuation and dissipation, the momentum distribution at zero 
temperature =< a*{k)a{k) > is 

oo 2^ e/^^ _^ {-2)Im[gMiik,uj + iO+)] 

-^^[zkSiu - Et) + (1 - zk)5{u: - E^)] = - 1, (24) 

where the label "Im" denotes the imaginary part. In deriving ()24p we have noted that 
quasiparticle dispersion E^ is always greater than or equal to zero and E^ is always smaller 
than or equal to zero. Because of this only the quasiholes give a contribution to the total 
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Figure 4: The momentum distribution in MI phases with fiUing factor uq = 2.0. Parameters 
are chosen as fohows: t/U = 0.01, fi/U = 1.5 and W/U = in (a) W/U = 2.0 in (b). 



density at zero temperature. Fig. 4 shows momentum distributions in MI phases and we 
find that three-body interaction W leads to not only flattening but also broadening of 
momentum distribution around k = 0. 

3.3 Superfluid Phase 

In terms of the superfluid phase, ipQ = {Gi^ + 4:t)/g. From the condition det{M2x2) = 0, 
we obtain four excitation branches ±Et^ with 



Compared to the Bogoliubov's theory for weakly interacting Bose gases, where only two 
excitation branches exist, there are two more excitation branches due to renormal- 

ization of the tunneling term by locally exact single-particle Green function Gic{iu}). Fig. 5 
show typical characters for E^: E'^ is gapped while i?^ is gapless. Moreover we also find 
that a gap opens between E^ and E^ when three-body interaction W becomes strong. 
Here, we are only interested in gapless excitations described by i-E^, which originate from 
spontaneously broken gauge symmetry and describe the low energy behaviors of the system. 
By expanding E^ in the vicinity of A; = 0, we find a linear spectrum i?^ = vk with 



.±2 



Ck 
hi 
h2 




(25) 



.1/2 



/i2 - hihs 



+ /l2 - /11/13 



V = 



a2 + 27(Gj+4t). 



Q = 



hi 



,7 = 



hi 



(26) 
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Figure 5: The excitation spectrum in superfluid phase with parameters t/U = 0.06, n/U 
1.5, W/U = 0.0 in (a) and W/U = 2.0 in (b). 



The existence of such hnear spectrum is consistent with the Landau's criterion of superflu- 
idity [24j . Fig. 6 shows the dependence of sound velocity v on three-body interaction W. 
It is should be noted that the system is in the superfluid phase near uq = 2 MI phase for 
parameters /U and t chosen in Fig. 6. For different parameters fi and t, on the one hand 
V generally is not monotonic as the function of VF, on the other hand the critical values 
Wc, where v = or the system enters into the MI phase, are very different. The second 
character can be illustrated from the phase diagram Fig.l: with the increase of W, the left 
boundary of MI lobe with hq = 2 changes very slowly, but the right boundary changes very 
quickly. 

From the excitation spectrum i-E^, the Green function gsF{k,iun) can be easily ob- 
tained 



gsFik,iuJn) 



iujn + 



+ 



lUJn 



where 



02 



+2 



aiEt =F ao 



+2 

k 



El 



El 



+ 



El 



+ 



Q- 



El 



(27) 



Ef^a^Ef 



aiE^ =F ao 



Ef) 



hi - hs + Ek + 2gVo'«i = ^3^i - /i2,ao = hi[h2 - /ii(efc + 29-00)]- (28) 



Following the same steps in (j24p . the momentum distribution in the superfluid phase is 

nk = Nijl6k,o-P+-Q+, (29) 

where the first term A^V'o^fc.o comes from the condensate and it is this term that supplies a 
coherent peak in absorptive images of atom gases. Fig. 7 show the momentum distribution 
from the second and third terms. In contrast to the situation in MI phases, three-body 
interaction W only broadens the width of the peak near k = 0. 



10 



(a) 



-t/U=0.07 
t/U=0.04 
■t/U=0.03 




-n/U=1.2 
- (i/U=1.5 
-^;U=1.8 



W/U 



W/U 



W/U 



Figure 6: The sound velocity v as the function of three-body interaction W for (a) ^/U 
1.2, (b) fi/U = 1.8, (c) t/U = 0.06. 





Figure 7: The momentum distribution in superfluid phase. Parameters are chosen as fohows: 
t/U = 0.06, ii/U = 1.5 and W/U = in (a) W/U = 2.0 in (b). 
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4 Conclusions 

In conclusion, an effective action for Bose-Hubbard model with two- and three-body on- 
site interaction in a square optical lattice has been derived by performing two successive 
Hubbard-Stratonovich transformations of the intersite hopping term. The main advantage 
of this method is that supcrfluid and MI phases can be analyzed on the same footing. 
Starting from this effective action, superfluid-MI phase transition, excitation spectrum and 
momentum distribution for two phases are calculated by taking into account Gaussian 
fluctuation about the saddle-point approximation. In addition the effects of three-body 
interaction are also discussed. We find that the sound velocity in superfluid phase generally 
is not monotonic as the function of three-body interaction, depending on the chemical 
potential and hopping term. 
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